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Abstract. We present deep Sofl and ISAAC near-infrared imaging data of the X-ray luminous galaxy cluster RDCS J 1252.9- 
2927. The ISAAC data were taken at the ESO Very Large Telescope under very good seeing conditions and reach limiting 
Vega magnitudes of 25.6 and 24.1 in the J- and X" s -bands respectively. The image quality is O'.^S in both passbands. We use 
these data to construct a colour-magnitude (C-M) diagram of galaxies that are within 20" of the cluster center and brighter 
than A" s = 24, which is five magnitudes fainter than the apparent magnitude of a L* galaxy in this cluster. The C-M relation 
is clearly identified as an over-density of galaxies with colours near J - K s = 1.85. The slope of the relation is -0.05 ± 0.02 
and the intrinsic scatter is 0.06 magnitudes with a 90% confidence interval that extends from 0.04 to 0.09 magnitudes. Both the 
slope and the scatter are consistent with the values measured for clusters at lower redshifts. These quantities have not evolved 
from z = to z = 1.24. However, significant evolution in the mean J - A" s colour is detected. On average, the galaxies in 
RDCS J1252. 9-2927 are 0.25 magnitudes bluer than early-type galaxies in the Coma cluster. Using instantaneous single-burst 
solar-metallicity models, the average age of galaxies in the center of RDCS J1252.9-2927 is 2.7 Gyrs. 
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1. Introduction 

The tight relation between colour and apparent magnitude for 
early-type galaxies in massive galaxy clusters (the C-M rela- 
tion) is seen at all redshifts, from the nearest clusters (Bower, 
Lucy & Ellis 1992 1 to the most distant clusters currently known 
(Rosati et al.[T333| Nakata et al. l2"0TJTl van Dokkum et al. l2Wl 
Stanford et al. l2"0T)2l Blakeslee et al. l2"0T)3l hereafter BFP). 

Observations show that although the zero-point of the C-M 
relation evolves considerably with increasing redshift (Aragdn- 
Salamanca et al. [T9^1 Stanford et al. nWgl hereafter SED; van 
Dokkum et al. EDTJTI Stanford et al. l^OTHl BFP), the slope of the 
relation and the scatter about it appear to evolve very little (Ellis 
et al. HW71 SED, van Dokkum et al. 1277001 BFP). However, 
for some clusters at z ~ 1, there is tentative evidence for a 
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flattening in the slope (van Dokkum et al. 2001 ; Stanford et al. 

i5om 

Complimentary studies of clusters up to z ~ 1 show that 
significant evolution is also occurring in the galaxy luminos- 
ity function (De Propris et al. Q998 Massarotti et al. 120031 
Toft, Soucail & Hjorth 2003 1 and the fundamental plane (van 
Dokkum et al. fT^98l van Dokkum et al. 120031 . Early-type 
galaxies in rich galaxy clusters are uniformly becoming both 
brighter and bluer as they become younger. 

The monolithic collapse scenario of Eggen, Lynden-Bell & 
Sandage (1962) is an attractive framework to model the ob- 
servations. In this scenario, the bulk of the stars form in a 
single burst over a relatively short period of time at redshifts 
greater than two (Ellis et al. 119971 Bower, Lucy and Ellis, 
1992 1. Hence, the scatter about the C-M relation and the slow 
steady evolution in the colours and luminosities of early-type 
galaxies are explained by the great age of the bulk of the stars. 
Extensions to this model allow future merging and star for- 
mation, but the scatter in the C-M relation limits the amount 
of merging and star formation that can take place (Bower, 
Kodama & Terlevich H998l . 
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This picture of passively-evolving, very old early-type 
galaxies in massive clusters should be compared to observa- 
tions of moderately distant clusters (z = 0.2 to z — 0.8) 
which show that the fraction of late-type, star-forming galax- 
ies in massive clusters increases with redshift (Dressier et al. 
P1W71 Couch et al. rrg38l van Dokkum et al. 120001 Nakata et al. 
2001 1, while the fraction of early-type galaxies does the oppo- 
site (Treu et al. 120031 van Dokkum et al. 2000 1. This picture 
should also be compared to semi-analytic and direct N-body 
numerical simulations which show that galaxy formation and 
evolution involves ubiquitous merging at all epochs (Kauffman 
and Chariot. fT^Sl Cole et al. 120001 Pearce et al. EOOTl . The 
simulations are able to reproduce the slope in the C-M relation 
and the scatter about it in present day clusters, although some 
difficulties, particularly at the high mass end, remain (Cole et 
al. 2000 1. At higher redshifts, the slope is predicted to flatten 
and the scatter is predicted to stay approximately constant, al- 
though there is a slight increase in the scatter at the bright- 
end of the C-M relation for clusters with z > 1 (Kauffman & 
Chariot. 1T9981 Ferreras & Silk l2000l 

Thus, there are two quite distinct pictures for the formation 
of early-type galaxies. In hierarchical merger models, the bulk 
of the stars form in disk-like galaxies that later merge to be- 
come early-type galaxies. In monolithic collapse models, the 
bulk of the stars form in early-type galaxies and subsequent 
merging and star formation are limited. In both models, the 
C-M relation is fundamentally a relation between the mass of 
a galaxy and the average metallicity of the stellar population 
(Faber [T973> . In hierarchical merger models, large ellipticals 
are formed from large spirals, which are better able to retain the 
metals that result from stellar evolution (Kauffman and Chariot 
1998 1. Similarly, in monolithic collapse models, larger ellipti- 
cals are better able to retain their metals. Although age differ- 
ences can be used to explain the slope of the C-M relation at 
low redshifts, the slope and the C-M relation itself are lost by 
z = 0.2 (Kodama & Arimoto fT997t if age is the sole reason for 
the slope. 

The morphological evolution that is seen in hierarchical 
models can lead to a bias (the progenitor bias) in morphologi- 
cally selected samples (van Dokkum et al. 2000 van Dokkum 
et al. 1200 1 i . The bias causes the progenitors of the youngest 
low-redshift ellipticals to drop out of morphologically selected 
high-redshift samples. Consequently, the C-M relation is simi- 
lar to that of a single-age stellar population formed at very high 
redshift and the scatter in the relation is approximately redshift 
independent. The progenitor bias is implicitly included in the 
semi-analytical simulations described above and allows the star 
formation history of early-type galaxies in clusters to be con- 
siderably more varied than that in monolithic collapse models. 
The model predicts that the fraction of early type galaxies in 
clusters decreases with increasing redshift. 

The importance of progenitor bias depends on the origin of 
the scatter in the C-M relation. If the scatter is entirely caused 
by age differences, then progenitor bias is important. If the scat- 
ter is partially caused by other effects, such as dissipationless 
merging with little subsequent star formation (van Dokkum & 
Ellis, 120031 or metallicity, then progenitor bias becomes less 
important and both the average age of the galaxies and the de- 



gree to which galaxies form coevally increase. The importance 
of progenitor bias also depends on the method used to derive 
C-M relations. C-M relations that are derived from morpholog- 
ical catalogues are more likely to be biased than C-M relations 
that are derived from photometric or complete spectroscopic 
catalogues. 

Although hierarchical models have become the standard 
model for describing the formation of early-type galaxies in 
both cluster and field environments, these models are unable 
to describe all the observational data. Whereas the hierarchical 
merging model predicts a dramatic difference in the star for- 
mation histories of early-type galaxies in the field and in clus- 
ters (Diaferio et al. 2001 1, only small differences are inferred 
from observational data (Willis et al. 2002, Treu et al. 2002 
Treu 2003 van Dokkum & Ellis 2003 1. More stringent tests of 
hierarchical models will come from observations of field and 
cluster galaxies beyond z ~ 1 . 

In this paper we describe deep near-infrared (NIR) obser- 
vations of RDCS J1252. 9-2927, an X-ray luminous cluster of 
galaxies at z= 1.237 (Rosati et al. 2003 1. These observations al- 
low us to construct a NIR C-M diagram of one of the most dis- 
tant massive clusters known to an unprecedented depth and ac- 
curacy. Throughout this paper, we assume £2m = 0.3, Q.\ = 0.7 
and Hq = 70 km/s/Mpc. In this cosmology, 1' on the sky cor- 
responds to approximately 0.5 Mpc at z — 1.237. Unless spec- 
ified otherwise, all colours and magnitudes are on the 2MASS 
system. 

2. Observations 

2.1. Soft NIR imaging 

As part of an NIR program to confirm high redshift galaxy clus- 
ters in the ROSAT Distant Cluster Survey (Rosati et al.[l998 1, 
RDCS J1252.9-2927 was observed during the nights of 1998 
July 9th, 10th and 1 1th with SofI (Moorwood, Cuby & Lidman, 
1998 1 on the ESO-NTT at the Cerro La Silla Observatory. SofI 
is equipped with a Hawaii 1024x1024 HgCdTe array, which, 
in the large field imaging mode, results in a pixel scale of 0729 
per pixel and a field of view of 4'9. The observations were done 
in J and K s , which, at the redshift of the cluster, approximately 
correspond to the rest frame V- and z-bands. 

Individual exposures lasted 10 seconds in K s and 20 sec- 
onds in J and six of these were averaged to form a single im- 
age. Between images, the telescope was offset by 10" to 30" in 
a semi-random manner. 

The data were reduced in the standard way. From each im- 
age, the zero-level offset was removed, a flatfield correction 
was applied, and an estimate of the sky from other images in 
the sequence was subtracted. Images with the best image qual- 
ity, defined here as the Full Width at Half Maximum (FWHM) 
of stellar sources, were then registered and combined. A sum- 
mary of the data is given in Tabled The central part of the SofI 
K s band image is shown in Fig. [2 

The atmospheric conditions at the time the data were taken 
were very good. Zero Points (ZP) were derived by observ- 
ing standards from the photometric catalogue of Persson et 
al. ([1998 1. Several standards were observed during each night. 
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Table 1. A summary of the observations taken with Sofl. The 
detection limit is the 5cr detection threshold over an l'/4 diam- 
eter aperture, which is approximately twice the stellar FWHM. 



Filter 


Exposure 


Image Quality 


Detection Limit 




(seconds) 


(") 


(Vega magnitudes) 


K, 


5400 


0.68 


21.1 


J 


5040 


0.72 


22.5 



The scatter in the J and K s ZPs throughout the entire run were 
less than 0.02 magnitudes. The Sofl J and filters are a 
good match to those used in the LCO (Persson) system, so 
no colour corrections between the natural Sofl system and the 
LCO (Persson) system are made. However, we do transform 
the magnitudes and colours to the 2MASS system (Carpenter 
2001 1 1 . Since the 2MASS and LCO systems are very similar, 
the transformations are small. 

We used the SExtractor software (Bertin & Arnouts 1996) 
to detect objects, to do the photometry and to classify sources 
as either point like or extended. The colours are derived from 
the flux within fixed apertures of 6 pixels (1'.'73). A small cor- 
rection (0.02 magnitudes) is applied to the J band data to ac- 
count for the slightly poorer image quality. The total magni- 
tude was estimated using the "BEST" magnitude in SExtractor 
(Bertin & Arnouts 1996 1. Since many of the sources in the cen- 
ter of the cluster are blended, the total magnitude was more 
often than not the corrected isophotal magnitude. 

The flux of non-stellar sources are corrected for galactic 
extinction. Using E(B - V) = 0.075 (Bouwens et al. I20031 . 
the corrections for J and K s are 0.067 and 0.027 magnitudes 
respectively. Stellar sources are not corrected for galactic ex- 
tinction. 

About a dozen sources from the 2MASS point source cat- 
alogue are visible in the Sofl images. The photometry for ob- 
jects brighter than J ~ 12.5 and K ~ 12 in the Sofl images 
is generally poor because of detector non-linearity. However, 
for fainter objects, we can make a comparison between Sofl 
magnitudes and those in the 2MASS catalogue. The variance 
weighted differences between the magnitudes in the 2MASS 
catalogue and the Sofl derived magnitudes for 9 objects in com- 
mon is AJ = -0.01 and AK S = 0.02. 

The C-M diagram of objects within 20" of the cluster cen- 
ter is shown in figure |2 An over-density of galaxies with 
J - K s ~ 1.85 is clearly seen in this diagram. Over the en- 
tire Sofl image other populations can be identified. Stars have 
J - K s colours that vary between and 1 .0, with a well defined 
peak at J-K s = 0.8, which corresponds to early M-dwarfs, and 
a less well defined peak at J - K s ~ 0.5, which corresponds to 
mid-G-dwarfs. 

The co-incidence of X-ray emission with two relatively 
bright galaxies (K s ~ 17.5) and the distinctive sequence of 
galaxies with J - K s ~ 1.85 showed that RDCS J1252.9-2927 
was probably a distant, rich galaxy cluster worthy of a more 
detailed study. We therefore initiated a program to obtain deep 
NIR images with ISAAC (this paper), deep optical and com- 
prehensive spectroscopic observations (Rosati et al. in prepa- 

1 See also http://www.astro.caltech.edu/2mass 
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Fig. 1. The central part of the Sofl K s -band image. The 3, 5 and 
7cr ROSAT X-ray contours are overlaid (Rosati et al. 1998 1. 
North is up and East is to the left. The image is 90" on a side, 
which, for the adopted cosmology, is 0.75 Mpc at z — 1.237. 



ration; Demarco et al. in preparation) and deep X-ray observa- 
tions (Rosati et al. l2003> . 

2.2. ISAAC NIR imaging 

RDCS J1252.9-2927 was re-observed with ISAAC (Moorwood 
et al. [19991 on Antu (VLT-UT1) at the Cerro Paranal 
Observatory. The observations were done in service mode and 
the cluster was observed on 25 separate nights, starting on 2001 
March 12th and ending on 2001 July 9th. 

ISAAC, like Sofl, is equipped with a Hawaii 1024x1024 
HgCdTe array, which, in the S2 imaging mode, results in a pixel 
scale of 07 1484 per pixel and a field of view of 2'5. The cluster 
was imaged in 7 S and K s . The J s filter is slightly narrower and 
redder than typical J filters and the blue edge of the filter is de- 
fined by the filter and not the atmosphere. This results in more 
stable photometry at the expense of introducing a small colour 
term. 

We did not center the cluster within the ISAAC field of 
view. Instead we used a series of pointings which put the clus- 
ter into the four corners of the array. The resulting union of 
images covers a region which is slightly smaller than the re- 
gion covered by the Sofl data, but is significantly deeper and 
has significantly better image quality. 

Individual integrations lasted 15 seconds in K s and 30 sec- 
onds in / s , and four of these were averaged to form a single 
image. Between images, the telescope was offset by 10" to 30" 
in a semi-random manner, and typically 30 to 40 images were 
taken in this way in a single observing block. In total, approx- 
imately 6 hours in both K s and J s were spent at each pointing. 
Since the cluster center was visible in all pointings, approxi- 
mately 24 hours is spent on the cluster in both 7 S and K s . 



4 



C. Lidman et al.: Deep near-infrared imaging of RDCS J1252.9-2927 




o - 



14 16 IB 30 22 

K s (2MASS) 

Fig. 2. The C-M diagram of objects in the Sofl data. The large 
symbols are objects that are within 20" of the cluster center 
and the small symbols are objects that are beyond 60" of the 
cluster center. Objects that have been classified as stellar are 
plotted with open symbols and non-stellar objects are plotted 
with solid symbols. No classification is made for objects fainter 
than K s = 20. The solid line is a fit to the C-M relation that was 
determined from the ISAAC data. 

In addition to these deep images, two additional regions 
that flank the eastern and western edges of the mosaic were ob- 
served so that potential cluster galaxies could be identified for 
later spectroscopic follow-up (Demarco et al. in preparation). 
The total exposure time in each of the these flanking fields is 
18 minutes, so the depth in these images considerably shal- 
lower that the depth obtained in the center of the mosaic. The 
details are given in Table 12 

The data were reduced in a similar way to that used for the 
Sofl data, but with following additional steps. 

- The difference in the relative level between odd and even 
columns was removed. The relative difference is a function 
of the average count level and it evolves with time. 

- An electronic ghost, which is most easily seen when there 
are bright stars in the field of view, was removed. 

- Sky subtraction is done with object masking. We used the 
XDIMSUM package for this step. 

- Images were corrected for field distortion, which can 
amount to several pixels at the edges of the ISAAC field 
of view. 

- Individual images were scaled to a common ZP before be- 
ing combined. 

The atmospheric conditions at the time the data were taken 
were very good. All but two of the nights were photometric 
and the image quality on the raw images varied between 0725 
to 078, with a median around 0745. Zero points were derived 
from the observations of photometric standards from the cata- 
logue of Persson et al. (1998). The ISAAC J s -K s colours were 



transferred to the LCO (Persson) system using the transforma- 
tion 

(J - K S ) LC0 = 1.028 * (J s - K S ) ISAAC - 0.011. 

As with the Sofl data, the colours and magnitudes were then 
transferred to the 2MASS system (Carpenter 2001 ) and only 
non-stellar sources are corrected for galactic extinction. 

The ISAAC data were taken over a large number of nights 
and, consequently, the image quality over the entire mosaic is 
slightly variable. The dispersion in the image quality as mea- 
sured from bright stars over the entire mosaic is about 8% in K s 
and 5% in 7 S . The median ellipticity is 0.05 and 0.03 for the K s 
and J s images respectively. In the central part of the mosaic, the 
uniformity of the PSF is considerably better, since the central 
part of the mosaic is common to all images. 

We used the SExtractor software (Bertin & Arnouts 1996 1 
to detect objects, to do the photometry and to classify sources 
as either point like or extended. The J-K s colour for objects in 
the central part of the mosaic are derived from the flux within 
fixed apertures of 6 pixels (0789 diameter). Since the image 
quality in the central part of the 7 S and K s ISAAC mosaics are 
very similar, there are no aperture corrections to the J - K s 
colours. The photometric errors were calculated independently, 
since the error in a fixed aperture is larger than the error that one 
would derive by multiplying the noise in a single pixel by the 
square root of the number of pixels in the aperture (Labbe et al. 
2003 1. The total magnitude was estimated using the "BEST" 
magnitude in SExtractor (Bertin & Arnouts 1996 1. 

The independently calibrated Sofl images are used to check 
the photometric accuracy of the ISAAC data. The difference in 
the absolute calibrations are less than 0.02 magnitudes and the 
scatter in the J-K s colour of relatively bright objects {K s < 19) 
is 0.02 magnitudes over a region that is within 30" of the cluster 
center. The scatter increases to 0.03 magnitudes if all objects 
within the region covered in Fig. [3] are included. The scatter 
increases with increasing area because image quality becomes 
more variable as one moves away from the cluster center and all 
colour measurements are based on aperture magnitudes. Unless 
otherwise stated, we restrict the measurement of cluster proper- 
ties in the analysis that follows to those galaxies that lie within 
20" of the cluster center. This region has been selected as a 
compromise between maximising the number of cluster galax- 
ies with respect to the number of field galaxies and obtaining a 
sufficient number of cluster galaxies for analysis. 

We also cross checked the photometry of the brightest clus- 
ter galaxies. Although there is significant dispersion in the total 
magnitude, about 0.15 magnitudes, the dispersion in the J - K s 
colour is less than 0.02 magnitudes and the mean difference is 
less than 0.01 magnitudes. The comparison between the inde- 
pendently calibrated Sofl and ISAAC data and the comparison 
between objects in the Sofl data and in the 2MASS catalogue 
suggest that the systematic error in our J - K s colours is no 
larger than 0.03 magnitudes. 

The central 75" of the ISAAC K s band image is shown in 
Fig. |3] This images can be compared to the image taken with 
Sofl (Fig.^, and it demonstrates the increased depth and spa- 
tial resolution of the ISAAC data. 
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Table 2. A summary of the observations taken with ISAAC. The exposure times and detection limits of the mosaic refer to the 
central part of the mosaic, where the exposures are deepest. In other areas, the exposure time will vary from one quarter to one 
half of this depending on the overlap. The detection limit is defined as the 5cr detection threshold over an 079 diameter aperture. 



Filter 


Region 


Exposure 


Image Quality 


Detection Limit 






(seconds) 


(") 


(Vega magnitudes) 


* s 


Mosaic 


81990 


0.45 


24.1 




Mosaic 


86640 


0.45 


25.6 




Eastern 


1080 


0.40 


21.5 




Eastern 


1080 


0.43 


23.3 


K s 


Western 


1080 


0.32 


21.5 


J, 


Western 


1080 


0.40 


23.3 




Fig. 3. The central part of the ISAAC /f s -band image. The image is 75" on a side, which, in the adopted cosmology, corresponds 
to 0.67 Mpc at z = 1.237. The two galaxies near the center of the cluster are 178 apart. Spectrally confirmed cluster members 
and field galaxies (Demarco et al. in preparation) are marked with arrows and circles respectively. Two arc-like features are also 
marked. 
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We mark the spectrally confirmed cluster members 
(Demarco et al. in preparation) that lie fully within Fig.[5]with 
arrows. With the exception of one cluster galaxy, which has 
the morphology of an edge on spiral and [Oil] emission, the 
morphologies of the spectrally confirmed cluster members are 
consistent with the morphologies of early-type galaxies. None 
of these galaxies show detectable [Oil] emission (Demarco et 
al. in preparation). For comparison, three spectrally confirmed 
field galaxies that lie with 20" of the cluster center are circled. 
The colours and morphologies of these three galaxies are quite 
different to the colours and morphologies of the spectrally con- 
firmed cluster members. One appears to be slightly disturbed 
edge-on spiral and the other two have highly irregular mor- 
phologies. 

There appears to be two relatively red arc-like features 
about 20" from the cluster center. Their shape and distance 
from the cluster are suggestive of giant gravitational arcs. 
However, higher resolution imaging and spectroscopy would 
be needed before conclusions could be drawn. 

3. Results 

3.1. The colour magnitude diagram 

The ISAAC C-M diagram of objects within 20" of the cluster 
center is shown in Fig. |4] Objects are plotted as solid sym- 
bols if they were classified as extended or as a star if they were 
classified as a point source. Objects fainter than K s < 21 are 
not classified and are plotted as solid symbols. Unlike the C-M 
diagram that was derived from the Sofl data (Fig. [2J and for 
reasons of clarity, objects outside this radius are not plotted. 
If they were, one would see a similar peak in the colours of 
stellar-like objects around J - K s ~ 0.8. 

Spectrally confirmed cluster members and field galaxies are 
marked with circles and crosses respectively. The field galaxies 
are also marked in Fig. [3] and all of them are morphologically 
distinct from the spectrally confirmed cluster members. 

Galaxies that are within 20" of the cluster center lie on 
a well defined C-M relation. The line in Fig. |4] is a fit to the 
galaxies within the blue rectangle (i.e. galaxies with K s < 21, 
J-K s > 1.5 and J -K s < 2.1) and within 20" of the cluster cen- 
ter. Even thought the fit has been done for galaxies brighter than 
K s = 21, the upper envelope in the colour of galaxies that are 
as faint as K s = 24 is defined by this relation. We also indicate, 
with boxes in Fig.0 the location of four spectrally confirmed 
cluster galaxies that are more than 20" form the cluster center 
but within the area bounded by Fig. [5] These cluster galaxies 
also lie on the C-M relation. 

Recently, Labbe et al ( 2003 1 measured the galaxy number 
counts to very faint NIR magnitudes. We have used their results 
to estimate that, to K s - 21, six field galaxies are expected 
within 20" of the cluster center, which is considerably fewer 
than the 36 galaxies we do observe. We do not separate field 
and cluster galaxies when fitting the C-M relation, since the 
three spectrally confirmed field galaxies are already excluded 
by our colour cuts, i.e. they lie outside the blue rectangle in 
Fig. |4] and the remaining small number of unidentified field 
galaxies are unlikely to bias our fit. 



The fit to the C-M relation was determined by adjusting the 
amount of intrinsic scatter that had to be added to the data until 
the reduced x 2 was one. A floor to the error in the colour, as 
determined by comparing the colours of the brightest galaxies 
in the ISAAC and Sofl data, was set to 0.02 magnitudes. Only 
those objects within the blue rectangle in Fig. @] were used in 
the fit. Separate fits were done for regions of different radii cen- 
tered on the cluster. In all cases, the dominant source of scat- 
ter is the intrinsic scatter. For completeness, we also report the 
fit to the slope of the C-M relation without adding additional 
scatter. The fits, the observed scatter and the inferred intrinsic 
scatter are reported in Table [3] 

In measuring the intrinsic scatter, it is imperative that mea- 
surement errors are accurately estimated. We checked the ac- 
curacy of our measurement errors by adding artificial galaxies 
directly to the reduced data and by processing the artificial and 
real data in an identical manner. The errors estimated in this 
way are in excellent agreement with the errors that are esti- 
mated from the image noise. Representative error bars are plot- 
ted in lower part of Fig. |4] for galaxies with apparent magni- 
tudes of K s = 18, 19, 20 and 21. 

We also used this technique to check for biases in the pho- 
tometry. There is a tendency for the SExtractor BEST magni- 
tude to miss an increasing fraction of the flux as objects be- 
come fainter. However, the error is less than 0.2 magnitudes 
over the magnitude range in which the fit to the C-M relation 
is done. Aperture magnitudes and the colours that are derived 
from them are unbiased. 

The slope of the C-M relation is -0.05 magnitudes per mag- 
nitude, which is quite similar to the slope measured in clusters 
with redshifts up to z ~ 0.9 (SED). If the fit is done using 
aperture magnitudes instead of the total magnitudes, the slope 
steepens slightly, although the change is within the measure- 
ment error. The scatter about the relation is unchanged. If the 
fit was repeated just for the spectrally confirmed cluster galax- 
ies that are marked in figure|3] the fitted slope and the measured 
scatter are also similar. 

Recent observations of clusters above z ~ 1 have pointed 
to a possible flattening in the slope (van Dokkum et al. 120011 
Stanford et al. 12002b . We see no evidence for a flattening in 
the slope of the C-M relation in RDCS J1252.9-2927. The dot- 
ted line in Fig.|4]is how the C-M relation of the Coma cluster 
would appear if it were moved to the redshift of RDCS J 1252.9- 
2927 (Rosati et al. 1999). Within the measurement uncertain- 
ties the slope for RDCS J1252.9-2927 is the same, but the J-K s 
colours of galaxies in RDCS J1252. 9-2927 are on average 0.25 
magnitudes bluer. This result is consistent with the trend seen 
in SED. 

We estimate the intrinsic scatter to be 0.06 magnitudes with 
a 90% confidence interval from 0.04 to 0.09 magnitudes. This 
is similar to the scatter seen in clusters with redshifts up to 
Z ~ 0.9 (SED). Even, if we have grossly overestimated the mea- 
surement errors, which is unlikely, the intrinsic scatter cannot 
be much larger than the measured scatter (0.06-0.08 magni- 
tudes). 

Within 20" (0. 16 Mpc) of the cluster center, 90% of galax- 
ies brighter than K s = 21, which is approximately 2.5 mag- 
nitudes fainter than the apparent magnitude of an L* galaxy 
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Fig. 4. The C-M diagram of objects within 20" of the cluster center. The diagram is generated from the ISAAC data. The solid 
circles represent objects that have been classified as extended and the stellar symbol represents the object that was classified as a 
star. Spectrally confirmed cluster members are circled and spectrally confirmed field galaxies are marked with an "X". The solid 
red line is a fit to the C-M relation of galaxies within the blue rectangle (see text) and the dotted line is where the E/SO sequence 
of the Coma cluster would lie if it were placed at z=1.24 (Rosati et al. I1999> . The dot-dashed lines are monolithic collapse 
models (Kodama & Arimoto 1997). From top to bottom they represent formation redshifts of z = 5, 3 and2. The error bars that 
are located near the bottom of the plot indicate the size of the colour and magnitude errors. Also plotted, as open squares, are the 
spectrally confirmed cluster members that are more than 20" from the cluster center but within the field of view of Fig. [3] 



in this cluster (Toft et al. in preparation) and 3.5 magnitudes 
fainter than the brightest cluster members, lie on the C-M rela- 
tion. Three of the seven galaxies that do not lie on the sequence 
were found to be non-cluster members from follow-up spec- 
troscopy. Of the other four, one is redder than the C-M relation 
and the other three are fainter than K s = 20. There is no pro- 
genitor bias in the central regions of this cluster. 



The average colour of a L* galaxy in RDCS J1252. 9-2927 
is J - K s ~ 1.85. This is similar to the colours measured 
for the early-type galaxies in RX J0848.9+4452 at z = 1.26 
(Rosati et al.[T999l, which have J - K ~ 1.85 on the UKIRT 
system. However, the transformation between the UKIRT and 
2MASS systems is uncertain, so it is not clear how compa- 
rable the average colour of galaxies in these two clusters are. 
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Table 3. The fit to the C-M relation. Only objects within the blue rectangle of Fig.0]are used in the fit. 



Radius 


Number of 


Slope 1 


Slope- 




Observed 


Intrinsic 


90% confidence region 


(") 


Objects 








Scatter 


Scatter 


for intrinsic scatter 


10 


14 


-0.046 


-0.058 


± 0.017 


0.066 


0.060 


0.039 to 0.091 


15 


19 


-0.052 


-0.055 


± 0.015 


0.068 


0.058 


0.041 to 0.086 


20 


29 


-0.046 


-0.048 


±0.015 


0.082 


0.070 


0.052 to 0.093 



Notes: 

1 No added scatter 

2 With additional scatter 



Using the transformation equations in Hawarden et al. (|7001 1, 
the colours agree very well, but using the transformation equa- 
tions of Carpenter (2001 1, one would find that the early-type 
galaxies RX J0848.9+4452 are 0.1 magnitudes redder. 

3.2. The epoch of galaxy formation 

We have used the monolithic collapse models of Kodama and 
Arimoto ( 1997 1 to illustrate how the C-M relation changes with 
formation redshift. In these models, star formation occurs dur- 
ing gas infall and terminates when the energy from supernovae 
ejecta exceeds the binding energy of the gas. Larger galaxies 
are able to hold onto their gas longer, and, consequently, star 
formation occurs over a more extended period of time and they 
have higher mean metallicities. Three models (with Z{ = 2,3 
and 5, where Z{ is the redshift of formation) are plotted in Fig. 
@]as the dot-dashed lines. They reproduce the measured slope 
very well and would suggest that the bulk of the stars in this 
cluster formed between z — 2 and z = 3. 

We have also used the simple stellar population (SSP) syn- 
thesis models of Bruzual and Chariot (2003) to estimate the 
mean age and the age spread of galaxies in the center of RDCS 
J1252.9-2927. We use J - K s - 1.85, which is the average 
colour of a L* galaxy in the cluster, and a spread in colours 
from J- K s = 1 .79 to J-K s = 1 .9 1 to estimate the mean age and 
the age spread. If we assume that the spread in colours is en- 
tirely due to age difference, then, for instantaneous single-burst 
models with solar metallicity, a mean age of approximately 2.7 
Gyrs and an age spread from 2.2 to 3.2 Gyrs are derived. These 
correspond to a mean formation redshift of Zf = 2.8 and a 
spread in formation redshifts from zt = 2.4 to 3.6. The up- 
per limits are highly uncertain. An increase of 0.01 magnitudes 
in J - Ks, which is smaller than our estimate of the system- 
atic error, will move the upper limit beyond Zf = 5. The model 
colours in these simple models do not get much redder than 
J - K s = 1.93 for solar metallicities. 

If these models are allowed to evolve passively to z = 0.5, 
the corresponding scatter in J - evolves to less than 0.01 
magnitudes, which is considerably less than the scatter mea- 
sured in SED. However, we caution that there are differences 
in the way C-M relations in this paper and Stanford et al. have 
been built. The C-M relation in this paper is created from the 
central 0.33 Mpc of the cluster, whereas those of SED are cre- 
ated from larger areas. With additional optical and spectro- 
scopic data on RDCS J1252. 9-2927, one will be able to extend 



the area over which the C-M relation is derived so that a more 
robust comparison can be made. 

Since the observed scatter in J - K s does not change from 
z = 1.237 to the present day, and since the scatter in SSP and 
monolithic collapse models should decrease as time goes on, 
it seems unlikely that the observed scatter in RDCS J1252.9- 
2927 is solely due to age differences. Alternative reasons for the 
scatter are metallicity variations at constant luminosity and/or 
dissipationless merging with little subsequent star formation. 
Such mergers may have already been identified in both cluster 
and field environments (vanDokkum et al. 1999 van Dokkum 
et al. l200Tl van Dokkum et al. l2003l . 

If most of the scatter is due to metallicity differences, the 
mean age of formation does not change, but the range of red- 
shifts over which galaxies would have formed narrows. If most 
of the scatter is due to dissipationless merging, then, in the 
monolithic collapse models of Kodama and Arimoto (19971, 
the mean age increases and the formation redshift is pushed 
higher, because galaxies will move left and down in the C-M 
diagram after the merging has taken place. The directions that 
galaxies move depend on the relative masses and metallicities 
of the progenitors. 

The derived formation epoch is directly related to the as- 
sumed metallicity and the star formation history. If a lower 
metallicity was assumed, an earlier formation epoch would be 
derived. This is the well known age-metallicity degeneracy, 
and the NIR imaging data that are presented here cannot break 
the degeneracy. An advance in this area will require a com- 
parison between deep spectroscopic data and the most recent 
stellar population synthesis models (Bruzual & Chariot 2003 1. 
However, in the SSP models that have been used here, the 
galaxies in RDCS J1252.9-2927 cannot be much metal poorer 
than solar, otherwise the formation epoch is pushed beyond 
the big bang. Additionally, the metallicity of the intra-cluster 
medium in this cluster is ~ 0.4 times solar (Rosati et al. 2003 1. 

4. Conclusions 

We have obtained very deep, J- and Zf s -band images of the X- 
ray luminous galaxy cluster RDCS J1252.9-2927 at z = 1.237 
with ISAAC on the ESO VLT and with Sofl on the ESO NTT. 
The data enable us to construct aJ-K s versus K s C-M diagram 
to K s = 24, which is five magnitudes below L* (Toft et al. in 
preparation). 

Galaxies within 20" of the cluster center define a tight C- 
M relation. The slope of the relation is -0.05 magnitudes per 
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magnitude and is similar to the slope measured in clusters at 
lower redshifts (SED). This strengthens the hypothesis that the 
slope in the C-M relation is due to metallicity and not age. We 
see no evidence for a flattening in the slope as predicted in 
hierarchical models and tentatively observed in clusters atz ~ 1 
(van Dokkum et al. 2001 Stanford et al. l^UOll . 

More than 90% of the galaxies within 20" of the cluster 
center and brighter than K s = 21 lie on the C-M relation. There 
is no progenitor bias in the centre of this cluster. 

The intrinsic scatter in the J - K s colour of galaxies about 
the C-M relation in RDCS J1252.9-2927 is 0.06 magnitudes 
and is similar to the scatter measured in clusters from z — to 
Z ~ 0.9 (SED). Hence, the scatter has not evolved from z = 1 .24 
to the present day. This weakens the hypothesis that the scatter 
in the C-M relation is solely due to age. Dissipationless merg- 
ing and metallicity variations at constant luminosity could also 
contribute to the scatter. We also see no evidence for increased 
scatter in the colours of galaxies at the bright end of the C-M 
relation. 

Our results can be compared to those derived from high- 
resolution optical images of RDCS J1252. 9-2927 that were 
taken with the Advanced Camera for Surveys on the Hubble 
Space Telescope in the F775W and F850LP filters. BFP find a 
tight C-M relation in the £775 - z&so versus zsso C-M diagram, 
and they show that neither the slope of this C-M relation nor 
the scatter about it have evolved from z — to z — 1 .24. This 
concurs with the findings of this paper. 

Using instantaneous, single-burst, solar-metallicity models, 
the average age of the bulk of the stars in the center of the 
cluster is 2.7 Gyrs. This corresponds to a formation redshift of 
Zf = 2.8. If the scatter about the CM relation is due to age, 
most of the galaxies in the center of this cluster were formed 
between z = 2.4 and z - 3.6. 
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